Astronomy & Astrophysics manuscript no. maris'etal 


February 2, 2008 


(DOI: will be inserted by hand later) 





The effect of signal digitisation in CMB experiments 

M. Maris 1 , D. Maino 2 , C. Burigana 3 , A. Mennella 4 , M. Bersanelli 2 ' 4 , and F. Pasian 1 



m 

O 
O 

(N 

Oh- 
< 



> 
On 
00 
O 
^1" 
O 
CO 

O 

S3 

Pl,: 

6 

C3 



1 INAF/Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 1-34131, Trieste, Italy 

2 Dipartimento di Fisica, Universita degli Studi di Milano, Via Celoria 16, 1-20131, Milano, Italy 

3 IASF/CNR, Sezione di Bologna, Via Gobetti 101, 1-40129, Bologna, Italy 

4 IASF/CNR, Sezione di Milano, Via Bassini 15, 1-20131, Milano, Italy 

Submitted to A&A April 1, 2003 

Abstract. Signal digitisation may produce significant effects in balloon - borne or space CMB experiments, when 
the limited bandwidth for downlink of data requires for loss-less data compression. In fact, the data compressibilty 
depends on the quantisation step q applied on board by the instrument acquisition chain. In this paper we present 
a study of the impact of the quantization error in CMB experiments using, as a working case, simulated data from 
the Planck/LFI 30 and 100 GHz channels. At TOD level, the effect of the quantization can be approximated as 
a source of nearly normally distributed noise, with RMS ~ q/y/!2N s , with deviations from normality becoming 
relevant for a relatively small number of repeated measures N a ~20. At map level, the data quantization alters the 
noise distribution and the expectation of some higher order moments. We find a constant ratio, ~ 1 / {s/V2o / q) , 
between the RMS of the quantization noise and RMS of the instrumental noise, o over the map (~ 0.14 for 
cr/q ~ 2), while, for <j/q ~ 2, the bias on the expectation for higher order moments is comparable to their 
sampling variances. Finally, we find that the quantization introduces a power excess, Cf x , that, although related 
to the instrument and mission parameters, is weakly dependent on the multipole £ at middle and large £ and can 
be quite accurately subtracted. For cr/q ~ 2, the residual uncertainty, ACp x , implied by this subtraction is of only 
~l-2% of the RMS uncertainty, AC" oise , on Cf v reconstruction due to the noise power, C"°" e . Only for &S30 
the quantization removal is less accurate; in fact, the 1/f noise features, although efficiently removed, increase 
CT Ise , AC"°' se , Cf and then ACf; anyway, at low multipoles Cf* > AC? oise > ACf . This work is based on 
Planck LFI activities. 

Key words, methods: data analysis - statistical; cosmology: cosmic microwave background 
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1. Introduction 

Since the successful detection of CMB anisotropy at angu- 
lar scales 9>7° by the COBE-DMR experiment (Smoot 
et al. 1992; Bennet et al.1996; Gorski et al. 1996), a 
large number of detections on smaller angular scales 
with increased sensitivity by both ground-based and 
balloon-borne experiments have followed (e.g. Lasenby 
et al. 1998, De Bernardis & Masi 1998, Bersanelli & Maino 
& Mennella 2002 for reviews). Recent exciting results from 
BOOMERanG (De Bernardis et al. 2000), MAXIMA-1 
(Balbi et al. 2000), DASI (Pryke et al. 2001) and CBI ex- 
periments (Mason et al. 2001) have provided strong sup- 
port to the inflationary scenario for structure formation 
and a universe with fip ~ 1- However to fully extract 
the wealth of cosmological information encoded in the 
CMB angular power spectrum, extreme mapping sensi- 
tivity is required. Only space missions can provide the 
proper combination of full-sky view and low level of sys- 



tematic effects that is necessary for "precision" CMB mea- 
surements (Danese et al. 1996), as demonstrated by the 
excellent results from the NASA satellite MAP 1 (Bennett 
et al. 2003). The ESA satellite Planck 2 (Mandolesi 
et al. 1998; Puget et al. 1998) is designed to measure the 
CMB anisotropy with an accuracy mainly set by astro- 
physical limits on a wide frequency range (30-900 GHz) 
and with unprecedent angular resolution and sensitivity. 

While space experiments benefit from unique environ- 
mental conditions, numerous technical issues need to be 
addressed. One of them is the limited bandwidth for down- 
loading of data to the Earth which often imposes strong 
limits on the available telemetry rate. For this reason it 
is often necessary to use loss less compression methods 
(Maris et al. 2000a) in order to reduce the data rate. 
Because the efficiency compression algorithms depends on 
the quantization step applied to the observed signal by 
the instrument acquisition chain on board the satellite, it 
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is necessary to find the optimal trade-off between the re- 
quired compression rate and the quantisation step, which 
needs not to be too high to preserve the scientific quality 
of the measured data. 

If the quantization step, q, is significantly smaller than 
the signal root mean square (RMS), ctt, and no satura- 
tion of the quantization levels occurs then the distortion 
induced by the signal quantization can be well approxi- 
mated by an additive source of non-Gaussian white noise 
with variance ~ g 2 /12 (Kollar, 1994). Corrections to this 
simplified model are required when ctt/q~ 1 , which depend 
on the kind of the total signal in input to the quantization 
stage of the on-board acquisition chain (Kollar 1994). For 
example, in the case of a pure white noise signal they are 
exponentially damped by a factor exp[— 2(7^<7T/g , ) 2 ]■ On 
the contrary, in the case of a sinusoidal signal (e.g. like 
the CMB dipole), corrections may by quite large and no 
simple analytical expressions can be written. 

It is worth noting that, in general, the observing strat- 
egy implies repeated observations of the same region of the 
sky on different time scales with detectors operating at dif- 
ferent frequencies. Therefore, considering that a sample is 
the result of the average of N repeated observations, we 
have that for a normally distributed signal with a^/q 3> 1 
the RMS of the quantization noise, a q is: 

(1) 



V12N 



In CMB experiments, however, the approximations 
of normal distributed signal and low a/q are often not 
strictly applicable. In fact the astrophysical signal is in 
general not normally distributed, and it is often required 
to use large quantization steps (i.e. <7t/<7 ~ !)■ The instru- 
mental noise power spectrum usually shows a 1// Q spec- 
tral behaviour, with a ~ 1 — 2 (see Burigana et al. 1997a 
and Seiffert et al. 2002), while the number of samples is 
not uniformly distributed when projected onto a sky map. 
Therefore it is apparent that the validity of eq. Q needs 
to be assessed in the context of space CMB experiments 
in order to evaluate possible deviations from this simple 
formula. 

In this paper we evaluate the impact of the quantiza- 
tion error on CMB measurements from space and provide 
tools and methods to optimize the quantization step to 
reach the required compression rate with minimal impact 
on the scientific output. 

We show examples relative to the Planck/LFI 30 
and 100 GHz channels. However, the formalism, tools 
and methods developed here are general and applicable 
to CMB anisotropy experiment based on differential re- 
ceivers, as well as to other kinds of experiments with large 
quantization steps. 

Using Monte Carlo simulations we estimate the effect 
of software quantization on the data considering various 
quantization functions and for different values of the quan- 
tization step. We also discuss the impact of different exper- 
imental conditions like the total number of repeated mea- 
sures for sample, N s , the white noise RMS, a, the l/f a 



noise and the scanning strategy. The impact of quantiza- 
tion is evaluated at three different levels: on time ordered 
data, on sky maps and their statistical moments, and on 
the angular power spectrum. 

In Sect. [3 we describe the analytical model of signal 
quantization and discuss the adopted numerical approach. 
Sect.[3]describes our results in terms of effects on TOD and 
maps including their statistical properties, and on the an- 
gular power spectrum. In Sect. 01 we present the possibility 
to reduce the effect of quantization in the recovered power 
spectrum. Finally, we discuss our main results and draw 
the conclusions in section |3J To improve readability, some 
mathematical details arc reported in Appendix lAl while in 
Appendix El we briefly verify that the extra-quantitazion 
step introduced via hardware by the analog to digital con- 
verter (ADC) quantization represents a minor effect with 
respect to the subsequent software quantization exten- 
sively discussed in the paper. In Appendix El we report 
on statistical moments of the map adopted in the simula- 
tions and of several foregrounds, for comparison. 

2. Analytical model and numerical simulations 

2.1. Analytical model 

During a typical CMB experiment from space, the mea- 
sured signal is first digitized and compressed on board and 
subsequently reconstructed on ground. Furthermore the 
scanning strategy usually involves multiple measurements 
from the same pointing direction that are subsequently 
averaged during ground data analysis. Therefore in order 
to assess the impact of the quantization error we need to 
consider the complete process, from quantization, to re- 
construction and averaging. In the following of this paper 
this process will be referred to as "QRA process" . 

Let us neglect the data compression step (that is out 
of the scope of our paper) and consider first the on-board 
quantization. In general the instrument output will be 
constituted by a Time Ordered Data stream (TOD) rep- 
resenting the sky temperature anisotropy; in the case of 
differential receivers, for example, the TOD will contain 
consecutive samples of temperature differences (sky mi- 
nus reference signal) so that the i th sample, ATj, will be 
quantized as (Maris et al.|5000a): 



AT? = Q 



AT, - AT, 



<l 



(2) 



where Q[x] and AT* represent the quantization operator 
and a zero level. 

The second step is the data reconstruction, i.e. the in- 
verse of eq. @. This second step is necessary because both 
the quantization step q and the offset AT* will not nec- 
essarily be constant during the mission; this implies that 
eq. (J5J) needs to be inverted in order to be able to compare 
measurements obtained in different periods of observation. 
The data reconstruction step can be defined as: 



A7f = (AT/ 3 + 0,) • q + AT, 



(3) 
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where C ( is an offset that depends on the considered quan- 
tization operator, and is defined by: 



O x =E[AT J ]/< ? -E[AT i Q ] 



(4) 



where E[Jf] represents the expectation of the random vari- 
able X. 

The value of O t is dependent on the quantization op- 
erator and assumes the values of 0, +1/2, —1/2 for the 
floor(x), ceil(x) or round(a;) operators, respectively. In 
the case of the trunc(x) operator the offset depends on 
the statistics of the random variable X; if X > then 
trunc(a;) = floor(x), if X < 0, trunc(a;) = ceil(a;). 

In our analytical treatment we will neglect the pres- 
ence of an offset; its effect will be discussed in the next 
section. Other simplifying assumptions are that the sig- 
nal is quantized in equally spaced steps and that AT* is 
always > 0. With simple considerations we can estimate 
the effect of the QRA process on the reconstructed maps 
and power spectra. We consider the case of typical CMB 
anisotropy experiments consisting of collections of many 
TODs each one based on the averaging over N s repetitions 
of the observation of the same stripe in the sky 3 . 

Assuming that the QRA acts mainly as an added white 
noise source, the effect of quantization on the power spec- 
trum can be easily estimated. Let C wr i T the white-noise 

induced excess of power and C^jy^ the equivalent excess 
of power induced by the QRA process; then we have that 
C %N S = K *n* 2 /N s and CI*** = K^ g y(12N s ), where 
N s represents the number of times that each pixel in the 
sky is measured during a single stripe and the pixel de- 
pendent normalization constants K^ n and K^ la are deter- 
mined by the scanning strategy (and the processing of the 
data, assumed to be a linear process). Since both quan- 
tized and not-quantized signals are sampled and processed 
in the same way, in principle we expect K^ n = K^ ra . 
Therefore, for any N s we have: 



QRA 



(5) 



The total noise power spectrum after quantization can be 
well approximated by: 



C total (~vwn 
e,N — ^i.n 



1 



(6) 



For a/q ~ 2 (as in the case of the Planck baseline) the 
noise power excess will be increased by ~ 2%. A similar 
relation holds for the ratio between the variance of the 
quantization noise and the variance of the white noise in 
a map: 



a 



map, total 



map, wn 



1 



12- 



q 



(7) 



3 For example, Planck TODs consist of about 1.1 x 10 4 
scan circles (for about 14 months of observations) each of them 
observed 60 times. See also section 2.2. 



where er m ap,wn is the RMS white noise in the map intro- 
duced by the instrument. This formula holds both for the 
noise average variance on the map and for the noise vari- 
ance of a given pixel in the map. 

2.2. Numerical simulations 

The analytical model described in the previous section al- 
lows a simple, first-order estimation of the QRA effect. We 
have performed numerical simulations to verify and quan- 
tify the validity and the accuracy of the analytical model 
and to study effects (like the presence of the skysignal 
and 1//" noise) not accounted for in the analytical ap- 
proach. The simulations are representative of the Planck 
mission. In particular we consider here the 30 GHz and 
100 GHz channels of Planck/LFI assuming angular res- 
olutions of ~ 33' and 10', respectively. 

2.2.1. Simulation tools 

To evaluate the effect of QRA process through numeri- 
cal simulations we have a) simulated the data acquisition 
during the space mission (with and without signal QRA 
processing) producing TODs; b) produced reconstructed 
sky maps and power spectra from simulated TODs; c) 
compared the results when including or not the QRA pro- 
cess. 

The simulation of the data acquisition during the 
mission has been performed using the Planck Flight 
Simulator developed by Burigana et al. Ijl997b|l . In this 
code we have included numerical tools to simulate data 
acquisition, signal quantization and signal reconstruction 
to allow the parallel generation of a quantized data stream 
plus a data stream of reconstructed data and to perform 
the statistical evaluation of the quantization error at TOD 
level (Maris et al. l2000b|l . 

The Flight Simulator is a code that simulates the CMB 
measurement according to the Planck satellite scanning 
strategy. Planck is a spinner (with the telescope observ- 
ing axis nearly perpendicular to its spin axis) that will 
map the sky through a large set of near great circles at 
a spin rate of 1 r.p.m. The satellite will orbit around the 
Lagrangian point L2 of the Sun-Earth system maintain- 
ing the spin-axis always in the anti-Sun direction. The 
measurement strategy is to acquire data from the same 
sky circle for about one hour before repointing the spin 
axis by approximately 2.5', so that the same sky pixel is 
measured ~ 60 times in each scan circle. 

The main simplifying assumptions concerning the 
scanning strategy that we have considered in our simu- 
lations are: (i) no drift or variation in the spin rate (equal 
to 1 r.p.m. for all the simulation), (ii) instantaneous re- 
pointing, (iii) perfect pointing with 85° observing angle 
with respect to the spin axis. 

The only astrophysical signal considered in the simu- 
lation has been a map describing the CMB fluctuations, 
therefore neglecting the effect of the dipole and of the 
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galactic and extragalactic contributions. To the astrophys- 
ics! signal we have also added white noise representative 
of the 30 GHz and 100 GHz radiometers (with an RMS 
of 1.056 mK and 3.325 mK, respectively, corresponding to 
integration times of 0.03 s and 0.009 s). In some simula- 
tions we also added a 1// component to study the effect of 
correlated noise on quantization; to this aim we considered 
two cases, with a 1// component having a knee frequency 
of 0.1 and 0.5 Hz, respectively. The first is representative 
of a typical worst case expected in LFI radiometers (for 
which we require a knee frequency < 0.05 Hz) while the 
second has been studied to enhance any possible effect 
connected to 1// noise but is clearly non representative 
of expected performances. 

The code developed by Burigana et al. I|1997b|l and 
Maino et al. ( 199SJ) has been used to generate maps from 
TODs. With this code it is also possible to significantly re- 
duce correlated systematic effects as the 1/ f a noise stripes 
from the data before generating the final maps. Power 
spectra have been obtained from the generated maps us- 
ing the anaf ast code of the HEALPix package from Gorsky 
et al. 1998, adopted here to pixelize the sky. 

2.2.2. Numerical procedures 

Let us now consider the time ordered data stream pro- 
duced by the Flight Simulator as an ordered sequence of 
measured temperature differences ATj, as in eq. (2). If we 
now take into account that for each pointing direction in 
the sky, indicated by the subscript p, on a given stripe 
(i.e., a scan circle for Planck) there will be N s repeated 
measurements then we can reorder the TOD according to 
two indices, (s,p), the first one (s = 1, . . . , N s ) represent- 
ing the repeated pixel measurements and the second one 
representing the pointing direction. Therefore we will have 
that ATi = AT SiP (i <-> (s,p)). 

Now for a single sample (s,p) the quantization error 
introduced by a quantization step, q, can be defined as: 

AT s E p = AT*, - AT S , P , (8) 

where the superscript R indicates a sample that has been 
quantized and then reconstructed. In a similar way we can 
define the following quantities averaged over iV s multiple 
measurements for each direction, p, in the sky: 

AT P = ±-X s AT SrP , (9) 

AT* = ^£ s AT s R p , (10) 

AT^ = AT* - AT p = is s AT^ . (11) 

To characterize the quantization error for the LFI sig- 
nal at TOD level we study the behavior of the distribution 
moments of AT p , AT* and AT^. In particular we ana- 
lyze their expectation, RMS, skewness and kurtosis as a 



0.1 - 




N 

s 

Fig. 1. Kurtosis for the distribution of the random er- 
rors for the QRA process as a function of N s . Crosses 
represents the kurtosis for the original simulated samples 
without quantization and reconstruction, but just averag- 
ing them in groups of N s . Bars represents the distribution 
of values of kurtosis for the quantization error for the aver- 
ages of N s quantized and reconstructed simulated samples. 
Note that the realisation of simulated TOD is not changed 
when q is changed so that the variability is just due to the 
QRA process. Circles represent the averaged kurtosis for 
the points forming the bar. The solid line represents the 
theoretical expectation for the kurtosis after N s averages. 

function of the quantization step q, of the quantization 
function Q[-\ and of the number of samples N s . 

The time ordered data with and without quantization 
have then be used to produce maps and power spectra in 
order to study the quantization error at these other two 
levels. 

The simulations have been carried out considering 
three distinct quantization steps, i.e. q w 0.5, 1.0, 2.0 mK 
for the 30 GHz channel and q w 1.2, 2.4, 4.9 mK for 
the 100 GHz channel. These three values of q correspond, 
for each frequency channel, to a/q ~ 2, 1 and 1/2, re- 
spectively. The case a/q ~ 2 corresponds to the current 
Planck baseline; the case a/q ~ 1/2 is representative of 
a case in which the noise model for the quantization dis- 
tortion is expected to fail, while the intermediate case was 
chosen evaluate possible deviations of the dependence of 
the quantization error on a/q from a power law. 

3. Results 

In this section we discuss the effect of signal quantization 
at the three different levels of time ordered data, final 
maps and power spectra. The results presented in this 
section have been produced considering the quantization 
operator floor(a;); this assumption can be made without 
any loss of generality, as all the quantization operators can 
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be considered equivalent provided that the proper offset 
is applied in the reconstruction formula, eq. OH 



3.1. Effect of quantization at TOD level 

To study the effect of quantization at the level of the time- 
ordered data stream we focus on the statistical properties 
of the data-stream by evaluating the various moments of 
the statistical distribution of the quantization error. 




3.1.1. Mean value of the quantization error distribution 

The first statistical quantity that we have analyzed is the 
mean value of the quantization error to evaluate if the 
QRA process introduces any spurious bias in the data and 
how this varies with q. This is relevant in the context of 
space missions, in which possible variations in the noise 
properties of the detectors will require slight modifications 
in time of q in order to maintain constant the a / 1 q ratio 
required for optimal compression. Therefore if any bias 
in the data is introduced by quantization it is reasonable 
to expect that this bias will vary in time (because of the 
changes in q) leading to systematic errors in the final maps 
and power spectra. 

We have calculated the mean value of the quantiza- 
tion error (AT 1 ! = mean[AT^ p ], see eq. for a ~ 1 
year-long data stream, considering different values of the 
quantization step q. Our results show that quantization 
introduces a bias in the data at a level jATgH/g < 10~ 3 
which implies |AT^| < 2 /j,K even for the largest values 
of q. This offset effect is small, and in addition it can be 
efficiently removed by destriping algorithms. 

Therefore we can expect systematic effects in the time 
ordered data at a very low level and the time scales of 
the updating of q on board (> 24 hours). Furthermore 
these effects can be efficiently reduced by applying de- 
striping and/or map-making methods to the time-ordered 
data and therefore they do not represent a concern 4 . 



3.1.2. RMS of the quantization error distribution 

The rms of the quatization error can be easily evaluated 
from the theoretical eq. . The validity of this relation- 
ship has been checked numerically by calculating the rms 
of the quantization error (indicated by Sn) for several 
TODs generated with various values of q and N around 
the reference values iVo = 60 and qo = 0.305 mK/adu, 
equivalent to a/q — 3.462 for the 30 GHz channel and 
a/q = 10.902 for the 100 GHz channel. A fit of S N with 
the relation: 



In Sn = In E + m N ln(N s /N Q ) + m q ]n(q/q ) 



(12) 



4 In fact, the current destriping codes are able to significantly 
reduce the impact of much larger drifts, such as those induced 
by the 1/ f°" noise or by periodic fluctuations induced by ther- 
mal fluctuations. 




2.00c-02 mK 



■ it ('AT B B0 )/RMS[J(|'AT E eo )] 



Fig. 2. The upper frame is the map of the quantization er- 
ror after the QRA process, plus destriping and map mak- 
ing, for a single 30 GHz radiometer, with a — 1.053 mK, 
a/q w 2, /knee = 0.1 Hz, iV s = 60. The lower frame is 
the ratio of the histogram of the map divided by the his- 
togram of a random variable normally distributed. The 
zero line represents the case of a quantization error nor- 
mally distributed. The y-scale represents the logj of the 
ratio. As indicated even in Tab. ^ the quantization error 
is not completely normal distributed. However the largest 
deviations are for large quantization errors, in excess of 
twice the RMS of the QRA process. 



yielded So ~ 9o/vl2iVo, ttin ~ — 1/2 and m q ~ 1 with 
an accuracy up to the third decimal figure. This confirms 
that in our case eq. JQ) is a very good approximation of 
the rms value of the quantization error distribution. 



3.1.3. Higher moments of the quantization error 
distribution 

In the ideal case of a normally distributed random signal 
with a/q 3> 1 the QRA error will also be nearly normally 

are 



5 In what follows we will use sometimes the terms "skew- 
ness" and "kurtosis" instead of skewness and kurtosis in- 
dex. These indices, combinations of the distribution statisti- 
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expected to vanish. On the other hand if a/q approaches 
1 then it is reasonable to expect a change in the values of 
the higher moments of the statistical distribution for the 
QRA error. 

From theory, the effect on higher moments may be 
evaluated from the characteristic function for the dis- 
tribution of the sampled values after QRA processing 
q-qQRA^-j _ $^( w )iJiv e (w), where ui is a real vari- 
able, <&Ar B (w) is the characteristic function for the QRA 
process (see appendix EJ, nAr B (w) is the correspond- 
ing function for the averaged input data. The expecta- 
tion for central moments of order k > being fik = 
{-i) k d k On%f A (u)/duj k \ ul =o- Since the QRA error is even 
distributed around zero expectations for its moments for k 
odd are null. If the same is true for the distribution of the 
unprocessed data then the expectations for all the /!& with 
k odd is null. The same does not hold for the expectation 
of moments with k even, in this case the QRA process 
may introduce a not null bias for some or all of them. In 
addition the QRA process randomly scatters about their 
expectation both odd and even moments. Then QRA pro- 
cess randomly changes the measured value of moments of 
any order k for a given realization respect to the values 
which have been measured in the case of no-QRA. Higher 
order moments, at least up to the fourth order, are needed 
to full-characterize the statistical properties of the QRA 
error. However the complete study of these features is be- 
yond the scope of this work and we limit our investigation 
to few cases illustrative of the methods which may be ap- 
plied to study them and the kind of perturbation which 
may be expected on real data. We may attempt a deeper 
investigation of this effect as numerical tools will improve 
enough to allow the generation and the analysis of a large 
data set of simulated missions. 

To evaluate the departure of the averaged quantization 
error from a normal distribution we have calculated the 
skewness and the kurtosis for the distribution of AT p , 
ATp and AT^ for various values of q and N s . 

As expected, the QRA process alters the sampling 
skewness in the input signal, so that it has meaning to 
put: 

ASkewness _ (Skewness[AT^] - Skewness [AT P ]) 
Skewness[AT p ] ~~ Skewness[AT p ] ^ 

In our simulation we observe values of 
|ASkewness/Skewness[AT p ]| up to 50% at the TOD 
level, comparable to the random fluctuations induced 
by the sampling variance. It is worth noting that this 
value depends on N s so that the variation in skewness 
decreases with increasing the number of samples of the 
same data point and that on the final map (in which each 
pixel will be the result of averaging a large number of 
measurement, especially close to the ecliptic poles) the 

cal moments [a, are defined by Skewness index = ^,3/ '/z^ 2 snd 
Kurtosis index = ^4/^2 — 3. 




sm(e)/sin(50') 

Fig. 3. Colatitudinal distribution of the variance of the 
QRA error on a sky map for a/q ss 2 and the 100 GHz 
channel. The straight line represents the best linear fit. 

impact on the skewness of the CMB distribution is even 
smaller (see Table |2J and varying with the colatitude. 

In Fig. n we show the behavior of the kurtosis for 
the QRA error (E [Kurtosis[ATp ] ] ) as a function of N s . 
From the plot it is apparent that the expectation of the 
kurtosis goes to zero for increasing values of N s . For an 
input signal with small expected kurtosis, the value of 
the sampling kurtosis at TOD level can be changed by 
the QRA process by a factor up to ~ 2, again compa- 
rable with the sampling variance for the kurtosis in our 
TODs. As for the skewness the impact on the kurtosis on 
maps is even smaller. However it is important to recall 
that, while skewness is changed randomly by the QLA, 
i.e. no bias is added by QRA to the skewness estimator, 
the same does not hold for the kurtosis. For the kurto- 
sis, the bias is -3/[l + 12(a/q) 2 } 2 , equivalent to ~ -0.001 
for a/q ~ 2. However, once the quantization step and the 
scanning strategy is defined, it is a simple matter to cal- 
culate this bias and to remove it. 

A non-zero value of the kurtosis in the distribution of 
the QRA error changes the statistical significance of the 
confidence limits respect to the usual definition of "stan- 
dard error" for which "lcr corresponds to a 68.27% confi- 
dence level" . To recover the standard error definition we 
have scaled the variance of the quantization error by a 
multiplicative factor apf B £l, which decreases monotoni- 
cally for increasing N B , whose lim.jv s — >oo &n s — 1> and 
ct\ = 1.18. Appendix 1X1 describes the computation of Q>n s . 

These results allow to safely apply the usual error 
propagation to combine both signal plus noise and the 
QRA error with an accuracy better than some /xK. The 
la confidence level for a given temperature average, /ijv B 
may be expressed as: 

^at s ± a Ns ± a Na (r q ,N e ~ Mjv 3 ± \J<j% b +a 2 Ns tf Ns . (13) 

The first term down the square root represents the effect 
of the white-noise, while the second term the effect of the 
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Fig. 4. Power Spectra from maps for simulated CMB, 
simulated CMB plus noise, simulated CMB plus noise plus 
QRA processing (labelled as Original, Not QRA, QRA) 
for the 30 GHz channel (1 radiometer). The inset is an 
enlargement of the power spectra for I > 400. For each 
power spectrum a smoothed spectrum is overlapped to the 
original spectrum. The a/q » 0.5 values has been chosen 
to increase the quantization effect. In the nominal case 
a/q w 2 the effect is about 16 times smaller and difficult 
to see in the plot. 



QRA error. Note that since the aw correction equations 
JHJ and J7J are valid only for N s sufficiently large to have 



\&n s — 1| "C 1. Otherwise the \y\2a/q) has to be scaled 
by a% in eq. © and eq. (0. In the case of Planck-LFI 
the measurement redundancy (N s ~ 60) is such that the 
distributions of the QRA error is sufficiently well approxi- 
mated by a normal distribution to apply the usual formula 
for the error propagation. 



3.2. Effect on the reconstructed map 

In Fig. |21 (top panel) we show a map of the QRA error 
for a/q « 2 for a 30 GHz detector. The map is pixelised 
according to the HEALPix scheme (Gorski et al. 1998) 
with a pixel-size of about 13.7 arcmin corresponding to 
A^sido (12 x A s 2 ide is the number of pixels in the map). To 
generate the map we first produced two TODs of the same 
simulation (including CMB, white and 1// noise) one with 
quantization and the other without. The error map was 
then obtained by differencing the maps generated by the 
two data streams. A simple "visual" inspection of the map 
indicates the absence of evident correlate structures 6 at 
few yitK level. Although the peak-to-peak variation of the 
quantization error on the map may appear large (~50 fiK) 
we must underline that quantization does not act as a 



6 The darker narrow band visible in the map is not due to 
quantization but to the overlapping of a subset of scan circles. 
In this region the integration time is greater and, therefore, the 
noise level is smaller. 



correlated systematic effect, but is mainly an added white 
noise at a few % level over the detector white noise (see 
eq. 0). 

In the bottom panel of Fig. |21 we show the ratio be- 
tween the QRA error distribution and a normal distribu- 
tion with the same mean and RMS, indicating that the 
QRA error is not completely normal distributed, with 
a number of error samples larger than 2a that is less 
compared to a normally distributed noise. These "non- 
Gaussian samples" , however, are mostly contained in a 
symmetric band at ecliptic latitudes between ±15°, and 
owing to the larger number of samples, this feature of the 
histogram is less large for the 100 GHz channel than the 
30 GHz. 

In Tab. ^ we report the % noise increase determined 
by quantization for various values of a/q. The comparison 
between the values determined by simulations with those 
calculated by eq. J7J shows the departure from the sim- 
ple theoretical model described in Sect. 12. ll For a/q ~ 2 
eq. (0) still represents a good estimate of the main effect 
of quantization. 

There is an interesting parallelism between what hap- 
pens to moments for a TOD of pure white noise when 
the QRA process is introduced, and what happens for the 
same moments in the map obtained from the same TOD. 
To illustrate it, let be to denote with fJ,k,map,wn,QR,A and 
l^k,TO D ,wn ,qra the fc-th moment for the map obtained 
from QRA processed data and its corresponding TOD 
from which the map is obtained, and with /ik,map,wn and 
Hk,TOD. wn the corresponding moments obtained for the 
map without quantization and the corresponding TOD. 
Looking to pixels as averages of N p samples, where N p > 1 
is the number of TOD samples entering a given pixel p in 
the map, it is possible to classify pixels as a function of 
N p . If [N p ] is the class of pixels averages of N p samples, 
g(Np) the partition function for the pixels over the various 
classes it is easy to demonstrate that, for any value of N s 
and a/q, and k — 2 the relation 



Mfc,map,iun,QRA _ A*fe, TOD, wn, QRA 



/^fc,map,'. 



Mfc,TOD, 



(14) 



holds exactly. From this eq. Q may be easily derived. For 
higher order moments (k > 2) the relation in eq. I|14l) holds 
better and better as N s is larger and larger. Indeed, if N s 
is not large enough other terms including g(N p ) will affect 
the right hand side of eq. Q. In the case of Planck hav- 
ing A s = 60 these terms are negligible at least up to k = 4. 
In conclusion, from the same reasoning leading to eq. JJJ 
it is possible to draw equations linking the variations of 
higher order moments on TODs to the corresponding vari- 
ations results of higher order moments on maps. 

As anticipated in the section related to the TODs, 
we compared the kurtosis in a map which has been 
QRA processed (Kurtosis [QRA]), against the kurtosis in 
the corresponding original map (Kurtosis [Orig]) giving 
the ratio of the variation (\AK/K\ = | Kurtosis [QRA] — 
Kurtosis[Orig]|/Kurtosis[Orig]). In general \AK\/K is 
greater for the 30 GHz channel than for the 100 GHz. In 
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Fig. 5. Power excess induced by signal quantization for the 30 GHz channel, 1 radiometer, i.e. the power spectrum for 
the quantization error map (panel a)) reported also in terms of STi — (Ce£(2£ + l)/47r) 1 / 2 (panel b)). Ratio between 
the power excess induced by the QRA process and the power spectrum of the unquantized map (noise plus sky) (panel 
c)) and its square root (panel d)). The approximate value of a/q for each curve is reported in panel a) (see Tab-^for 
the corresponding exact values). 



particular for the 30 GHz channel considering maps which 
does not contain 1 // noise and consequently have not been 
destriped the kurtosis may be varied even of a 30% by the 
QRA process. However, comparing maps which does not 
contain 1// noise and have not been destriped with maps 
originated from data containing 1/f noise and destriped 
it is possible to see that the relative effect on the kur- 
tosis is much smaller than in the previous case. This is 
nothing else than the fact that 1/f noise and destriping 
increases the kurtosis of the unquantized signal given at 
the denominator of the \AK/K\ ratio, but leaves prac- 
tically unchanged the AK difference. From Tab. Q it is 
also possible to see how \AK/K | scales with a/q both for 
the 30 GHz and 100 GHz maps. Comparing maps which 
contains even the 1/f noise and have been destriped (last 
column of \AK/K\M.2% (1.5%). For both 30 GHz and 
100 GHz maps \AK/K\ increases up to a factor of about 
7 decreasing a/q from 2 down to 1/2. 

The square root of the sampling variance of K is 
y/2A/N~ (see e.g. Kendal & Stuart 1977), i.e. V2/N sidc 



in the case of the HEALPix scheme, over a map with 
A^pix pixels obtained from a single receiver (radiome- 
ter for LFI). For maps at N s i dc = 256 (512) this gives 
5.5 x 10~ 3 (2.8 x 10~ 3 ), which is equivalent to a sensitiv- 
ity \AK/K\ ~ 30 - 50% (20%) over a map from a single 
receiver in the case of a sky of pure CMB fluctuations (see 
Tab. C.l in Appendix C). When \AK/K\ is scaled to take 
into account that the final map is derived from a combi- 
nation of more radiometers, the sensitivity to K improves 
to 0.0028 (0.00057) at 30 GHz (100 GHz), equivalent to 
\AK/K\ ~ 15-25% (4%). The \AK/K\ values induced by 
the QRA process (see Tab. 2), although not large, is not 
negligible in the case of maps of pure CMB anisotropy plus 
noise, particularly because it is not a statistical but a sys- 
tematic effect. At 30 GHz, the single receiver noise is small, 
and the above conclusion does not change significantly in- 
cluding also noise contribution to the kurtosis, while at 
100 GHz, where a single receiver noise is quite large for 
LFI, the reference value of K changes significantly includ- 
ing or not the noise, and the comparison between \ AK / K\ 



M. Maris et al.: Digitisation error in CMB Experiments 9 

30 GHz 

a/q % variance increase % variance increase a b 

(from simulations) (from eq. (mK 2 ) (mK 2 ) 

2.076 L94 L93 (-9.00 ± 1.16) x 10" Y (2.105 ± 0.011) x 10~ b 

1.038 7.67 7.73 (-4.02 ± 0.46) x 10~ 6 (8.459 ± 0.046) x 10 -5 

0.519 30.94 30.59 (-1.30 ± 0.19) x 10~ 6 (3.399 ± 0.019) x 10^ 4 



100 GHz 



a/q 


% variance increase 


% variance increase 




a 


b 




(from simulations) 


(from eq. |7J 




(mK 2 ) 


(mK 2 ) 


2.724 


1.13 


1.12 


(- 


-1.436 ± 0.036) x 10~ b 


(1.5061 ±0.0036) x 10 -4 


1.362 


4.55 


4.49 


(- 


-5.903 ±0.146) x 10" s 


(6.0464 ± 0.0145) x 10 -4 


0.681 


17.57 


17.97 


(- 


-2.241 ±0.056) x 10" 4 


(2.3698 ± 0.0055) x 10 -3 



Table 1. Quantization error at map level. The table compares the % variance increase determined by simulations 
((^map total /""map.wn — x ) x w ith the expectation derived from eq. 10. The last two column report the best fit 
parameters to estimate the RMS quantization error as a function of the pixel colatitude. 



30 GHz 



a/q 


1 


Skewness(X) 


Kurtosis(X) 


\AK/K\% 


\AK/K\% 




[mK] 






(no i/n 


(1/r + dart.) 


2.076 


0.509 


-0.00029069 


0.370293 


30.7 


4.2 


1.038 


1.017 


0.00692233 


0.394023 


62.8 


5.9 


0.519 


2.035 


0.00547806 


0.383053 


63.9 


26.1 








100 GHz 






a/q 


q 


Skewness(X) 


Kurtosis(X) 


\AK/K\% 


\AK/K\% 




(mK) 






(no 1/D 


(1/r + dent.) 


2.724 


1.221 


0.00333290 


0.257705 


0.4 


1.5 


1.362 


2.441 


0.00031871 


0.274681 


3.1 


5.4 


0.681 


4.883 


-0.00492754 


0.267852 


14.9 


8.8 



Table 2. Table of sampling statistics for the quantization error maps: skewness and kurtosis for the QRA. The fifth 
column lists values of \AK/K\ x 100 obtained comparing the kurtosis of maps not containing 1// Q noise and not 
destriped, with and without QRA processing. The same for the last column where l// a noise has been introduced 
and destriping has been applied. For other details see Tab. Q 



representing the sensitivity and that reprenting the QRA 
effect may depend strongly on the considered realization. 

On the other hand, the microwave sky includes 
fluctuations also from Galactic and extragalactic fore- 
grounds, with remarkable non Gaussian distributions. In 
Appendix C we report the statistical moments of the 
most relevant microwave components and of composite 
maps. As evident, even including Galactic cuts and re- 
moving bright sources and Sunyaev-Zeldovich effects (see 
last three lines of Tab. C.3), foregrounds change the 
sky kurtosis index with respect to the case of a pure 
CMB anisotropy sky much more than the QRA process. 
Therefore, the foreground removal is an aspect much more 
important than QRA process for the estimation of CMB 
anisotropy high order moments. 

Finally in Fig. [3] we show the variation with the ecliptic 
latitude of the variance of the QRA error, which decreases 
as a function of sin 8, where 6 is the pixel colatitude. The 
figure also shows that for 6 > 15° the variance can be well 
approximated with a + 6 sin (9/ sin 50°, where the best-fit 



coefficients a and b are listed in the last two columns of 
Tabled This relationship represents a useful method to 
estimate the rms QRA error at a given latitude. 

3.3. Effect on the power spectrum 

In Fig. 0] we show an example relative to a Planck-LFI 
30 GHz radiometer that highlights the qualitative fea- 
tures of the quantization effect on the measured power 
spectrum. To enhance the effect we have chosen a value 
a/q ~ 0.5, which is clearly not representative of typical 
conditions. 

The figure shows three power spectra. The first one 
(lower curve) is the power spectrum of the original map. 
The second (middle curve) is the power spectrum of the 
map "observed" by the simulator without quantization; 
the power increase at high multipoles is due to the instru- 
mental white noise. The third (upper curve) is the power 
spectrum of the "observed" map including quantization. 
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30 GHz 



a/q 


<7 


fiance 


prec. 


/ /1CXC \ 

\ W>70,60/ 




(mK) 


(Hz) 




{W- g rnK 2 ) 


2.076 


0.509 


0.1 


no 


0.331 ± 0.020 


1.038 


1.017 


0.1 


no 


1.323 ± 0.081 


0.519 


2.035 


0.1 


no 


5.410 ± 0.350 






100 GHz 




a/q 


q 


/knee 


prec. 


^QRA 




(mK) 


(Hz) 




(W- g rnK 2 ) 


2.724 


1.221 


0.1 


no 


0.561 ± 0.027 


1.362 


2.441 


0.1 


no 


2.242 ± 0.110 


0.681 


4.883 


0.1 


no 


8.877 ± 0.470 








30 GHz 




a/q 


q 


/knee 


prec. 


^QRA 




(mK) 


(Hz) 




{W- g rnK 2 ) 


1.038 


1.017 


0.5 


no 


1.323 ± 0.076 


0.519 


2.035 


0.1 


yes 


5.820 ± 0.350 



Table 3. Values of C^, A obtained from simulations for 
a single 30 GHz and 100 GHz horn. Errors represents 1 a 
dispersion of Ce about the average. The last subtable re- 
ports the effect of a change of /knee for the baseline scan- 
ning strategy and for a scanning strategy which includes 
a precession of the spin axis. 



The inset is an enlargement of the power spectrum for 
400 < £ < 767. 

From the inset it is evident that the QRA process in- 
troduces a power excess in the spectrum. A closer look 
reveals that most of the features of the original power 
spectrum are preserved. At multipoles larger than a crit- 
ical value 7 (£ > £ C iit) this power excess approaches a 
value, C^ RA , nearly independent on I and dependent on 

N s (we use here N s — 60). The value C^ RA can be used 
as a convenient estimator of the quantization effect on the 
power spectrum at large multipoles and can be quite well 
estimated by the average of the power excess for £ > £ CY it- 



C 



QRA 
W,60 



E 



max(<) ^yQRA 
it °A60 



max(£) - 4 r it + 1 



(15) 



where max(f) ~ 750 and 1500 respectively for the 30 and 
100 GHz frequency channels. 

The power spectra of absolute and relative quantiza- 
tion error maps, in terms of both Ci and STg are shown in 
Fig. [S] for different values of a/q. The power excess varies 
roughly as a power law of cr/q, in good agreement with 
eq. ©. 

Table 01 gives C^/^q for three quantization steps for a 
single radiometer at the 30 and 100 GHz frequency chan- 
nels. The simulations from which Cft 60 are derived as- 
suming the nominal PLANCK same scanning strategy 
and l/f noise knee frequency. In addition we report two 



7 4rit ~ 70 for the 100 GHz channel and i cilt 
30 GHz channel. 



30 for the 



cases obtained changing the knee frequency and the scan- 
ning strategy, which will be discussed at the end of this 
section. 

A fit of the data obtained by simulations with eq. iJBJ 
yielded a residual log- log correlation coefficient |1 — 
Hog-log | < 10~ 5 , which indicates an excellent agree- 
ment with the expected q 2 dependence. This allows to 
parametrize 8 : 



c: 



QRA 



c: 



QRA 1 



(16) 



The normalization factor can be easily deter- 

mined from the above-mentioned best-fit; a good approxi- 
mation can be obtained by replacing C^ RA with Kf / A P i X , 
where A P i X represents the number of pixels on the map 
and Kf is a normalization factor 0.72 < Kf < 1 that 
depends on details such as: the scanning strategy, the pix- 
elization scheme and the geometry of the instrument. A 
further approximation Kf = 1 yields: 



c: 



QRA 



UN S N„ 



(17) 



which is a useful relationship to estimate upper limits to 
the quantization effect on the recovered power spectrum. 
Finally we can determine the ratio K^, n /K 2 Ya defined 

in eq. © by fitting logC^ A versus log (Vl2a/q) 2 and 
taking the ratio with the power spectra on the unquan- 
tized noise. This yields K^/K^ = 1.024, that represents 
a good approximation of the ideal value K%, n /K* I& = 1. 
This is a further confirmation that eq. (J3J) and eq. © 
well describe the quantization effect on the power spec- 
trum. From these simulations the total power excess due 
to white noise plus the QRA processing is: 



wn; 



3.55- 



(18) 



On the other hand, as evident in Fig. 0] referring to 
realistic simulations including the l/f noise and the ap- 
plication of the destriping algorithm, for £ < £ CI it, the 
power excess is no longer constant, but shows an increase 
for decreasing I. However, for a/q « 2 the deviation from 
a constant power excess at low £ is not large. For exam- 
ple, for the 30 GHz channel it is equivalent to <16% in- 
crease of the power excess (50% for the 100 GHz). For 
a/q w 1 the deviation is significant up to 7.5C^ RA at 
30 GHz and 16.8C £ QRA at 100 GHz. Although It is not easy 
to parametrize this effect, C^ RA may be well described by 
a second or third order polynomial of log £ in the log-log 
space (Maris 120023 for £ > 1, which for £ > 70 well fits 
the constant power excess previously assumed. However 

8 In this work we have not studied numerically (because of 
computation time limits) the dependence of C^ RA on N s . So 
the 1/JVs dependence of eq. © Cf RA as a function of a, a/q, 
N B has been assumed. 
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the results of this fit can not be easily generalized as in the 
previous case, since details such as the polynomial degree 
are related not only to the mission parameters but also 
to the assumptions implicit into the destriping procedure 
and the parameters of the 1// noise. This is supported by 
the comparison of simulations obtained for different val- 
ues of /knee (see Tab.0]for the 30 GHz frequency channel). 
No significant differences appear for i > 30 by changing 
/knecfrom 0.1 Hz to 0.5 Hz, while at I < 30 the power 
excess changes up to a factor of two. 

A change in the scanning strategy instead, affects the 
power excess at both low and large I. In Tab. 01 is shown 
what happens allowing for a precession of the spin axis 
during the mission. C^ RA increases by about 7.6% at high 
£, and as much as 82% at low I. Since the same sequence of 
random noise samples has been used for both the precessed 
and un-precessed scanning strategy, these changes are due 
to the different weighting of the samples in the sky induced 
by the precession of the spin axis. 

4. Removing the effect of the QRA error in the 
power spectrum estimation 

We have seen that quantization effects on the power spec- 
trum are significantly smaller than the noise and that they 
are well represented by few parameters. Here we discuss 
the possibility to remove quantisation excess noise in the 
data analysis from the final estimated power spectrum. Of 
course, this is of particular interest at middle and large 
multipoles while at low multipoles the noise power spec- 
trum is much smaller. The estimation of the sky power 
spectrum can be derived from the power spectrum ob- 
tained from the observed map by subtracting the sum of 
the expected noise power spectrum and of the expected 
QRA power spectrum: 

Cf y = C,° bs - E[C; oise ] - E[C^ QRA ] . (19) 

The accuracy of the subtraction of the QRA effect is deter- 
mined by the accuracy of the model for C^ RA and by the 
dispersion of the true values of C^ RA about their expecta- 
tion. The typical accuracy of this subtraction at large £, 
as evaluated from the last column of Tab. [21 is of about 
5% (6%) for the 100 (30) GHz channel, leaving a residual 
(i.e. after subtraction) QRA error equivalent to an unsub- 
tracted QRA error with a quantization step at least four 
times smaller than that really used. 

The RMS of the expected QRA power spectrum de- 
creases with £, analogously to the case of the RMS of 
white noise that exibits a typical 1/v^ dependence. We 
expect then that the accuracy of the subtraction of the 
QRA effect improves with increasing I. For example, tests 
with the simulated power spectra reported in section [3] 
show that the RMS of the QRA power spectrum ranges 
between - lO" 11 mK 2 and - fewx lO" 10 mK 2 at I ~ 1500 
(750) for the 100 (30) GHz channel when a/q ranges be- 
tween 2 and 1/2. As a consequence, the power excess is 
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reduced up to a factor 45 (34) for the 100 (30) GHz chan- 
nel at very large £, equivalent to an unsubtracted QRA 
error with a quantization step reduced by a factor ~ 6 — 7. 
On the contrary, at I of few tens the RMS of the expected 
QRA power spectrum is of the same order of magnitude of 
the power spectrum itself, Cfgg. For example at I ~ 15 it 
ranges between some x 10 -10 mK 2 and some x 10~ 8 mK 2 
for the 100 GHz channel, and between few x 10 -10 mK 2 
and few x 10 -8 mK 2 for the 30 GHz channel. Compared 
with the averaged power excess for ^<30 this is equivalent 
to a reduction of a power excess between ~ 50% and a 
factor of two. 

We searched for possible improvements of the removal 
procedure by with a polynomial fit of log Cf^ versus log £. 
We find that the best fitting relation depends on the fre- 
quency channel, the value a/q, the scanning strategy and 
the beam position so that a proper fit for each receiver and 
scanning strategy would be required. However, our simu- 
lations show that for a/q m 2 the values of Cfgg are very 
dispersed around the mean for I < 30, while for I > 70 
the goodness of fit does not depend significantly on the 
degree of the polynomial. Therefore a simple linear fit is 
adequate. 

In order to evaluate the final impact of the accuracy 
of the subtraction of the QRA quantization effect in the 
final estimation of the sky power spectrum it is crucial to 
compare the RMS of the power spectrum of QRA quan- 
tization error with the RMS of the noise power spectrum 
that sets the fundamental uncertainty in the sky power 
spectrum recovery. 

We require that 

RMS [C$§] < RMS [Cf oise ]/A , (20) 

where the constant A has to be set large enough to reduce 
the residual QRA quantization error to an acceptable level 
(for example, A > 5 implies a residual QRA quantization 
error less than 20% of the unsubtractable instrumental 
noise) . 

After having applied the destriping algorithm, the in- 
strumental noise is essentially white noise dominated and 
we can then assume 

RMS[Q noiso ] ~ RMS[C ; WN ] ~ E[C ; WN ]/\/Z (21) 

As discussed before, Cf go depends on the considered 
receiver and scanning strategy; of course, the same holds 
for RMS[C ; c gQ]. It is then impossible to rewrite in general 
the condition (|20|l in terms of a/q. 

To circumvent this problem we have searched for a 
law, Ui, which gives at each I the maximum among the 
ratios RMS[Cj c gQ]/Cj c gQ found for the set of simulations 
considered here. We find that the approximate law 

Ut ~ OMl' 057 (22) 

is a quite good approximation for the whole range of mul- 
tipoles relevant here. 



12 M. Maris et al.: Digitisation error in CMB Experiments 

Frequency B A A" C £ Q f A K f 

Channel {\og 10 mK 2 ) (log 10 mK 2 ) (mK 2 ) 

30 GHz +0.015 ±0.022 -8.8425 ± 0.0013 -8.8903 ± 0.0156 (0.7724 ± 0.0120) x 10"' 0.7241 ± 0.0113 
100 GHz -0.008 ±0.023 -8.3834 ± 0.0018 -9.4257 ± 0.0125 (0.2251 ± 0.0028) x 10 -7 0.8364 ± 0.0105 

Table 4. Results of the fit of log 10 C^ Q A versus log 10 a/q for the data in Tab. [3] Only the data in the first and the 
second subframes have been included in the fit. 



The condition expressed by eq. (f^Uf) is then satisfied 
provided that 



0.93f 



-0.57 



E[cfi c }<E[cr N }/V~e/x. 



(23) 



At middle and large multipoles the ratio between the 
variances of the white noise and of the quantization error 
equals that between the two corresponding power spectra. 
By using eq. J7J), the condition 123(1 can be then easily 
rewritten in terms of the a/q ratio as: 



0.93 
'12-' 



-0.07! 



(24) 



At £^30 (500) this condition is clearly satisfied even 
for A>4, 16 or 65 (5, 20, 80) respectively for a/q ~ 1/2, 1 
or 2. 

In conclusion, this analysis demonstrates that it is pos- 
sible to subtract the quantization error impact in the sky 
power spectrum estimation from Planck data at a level 
of accuracy which makes its residual effect very small and 
significantly smaller that the intrinsic unsubtractable un- 
certainty due to the noise. 

However it should be noted that these results are ob- 
tained with the assumption of stationary noise. In partic- 
ular the application of this correction requires to consider 
every changes in the detector calibration which affects q 
along the mission. Moreover, the previous result about 
the sensitivity of C^ RA on the scanning strategy imposes 
to consider that the final C^ RA to be subtracted from 
the power spectrum shall be evaluated by numerical sim- 
ulations to take into account the real scanning strategy. 
This is particularly relevant for a multi-horn instrument 
like Planck, where a full evaluation of C^ RA shall be 
obtained only with a full simulation of the mission with 
each horn observing somewhat differently the sky. 

5. Conclusions 

An unprecedented amount of high quality data are ex- 
pected from the new generation of CMB experiments de- 
signed for extended, high resolution imaging of the mi- 
crowave sky. These data sets will be analyzed to measure 
the angular power spectrum to high precision, and also 
to search for low amplitude statistical properties of the 
CMB distribution, such as non-gaussianity or high-order 
autocorrelation functions. It is therefore of interest to ac- 
curately consider the effect of instrumental systematics 
that could significantly perturb the statistical properties 



of the noise and signal distribution both in the map and 
in the angular power spectrum recovery. 

For a space mission, the need of minimizing the teleme- 
try rate requires that the signal digitization in terms of 
the a/q ratio be as low as possible. On the other hand, 
a poor signal sampling, combined with the adopted scan- 
ning strategy, may produce subtle effects impacting the 
statistics. In this paper we have analyzed the impact of 
signal digitization, reconstruction and averaging of a full- 
sky CMB survey. The effect of the QRA process has been 
analyzed at the level of data streams, maps and power 
spectra, both by analytical approximations and numerical 
simulations. We have simulated the effect for the case of 
the Planck/LFI experiment, but the methods presented 
here and the main conclusions can be applied to any other 
experiment. The simulated signal includes white and 1//- 
type noise and a realization of CMB anisotropy, while the 
effect of QRA process on high order statistics has been 
compared with that from main Galactic and extragalactic 
foregrounds. 

Because the signal is dominated by white noise, the 
a/q ratio is the main parameter describing the QRA 
scheme. At the level of data streams we studied the dis- 
tribution of QRA errors as a function of a/q and of the 
number N s = 3, 4, . . ., 60 of repeated pixel observations 
obtained by the PLANCK scanning strategy. As expected, 
when N s is large the QRA error is well approximated by a 
normally distributed random variable, and standard error 
propagation may be applied. 

For large N s most of the numerical results at the 
data-stream level may be well approximated by analyt- 
ical formulae. At map and Cg level this is not always true, 
since the effect of the scanning strategy combined with 
destriping and map-making algorithms can not be mod- 
elled in a simple way. For these cases, we have derived 
semi-analytical relations whose parameters are calibrated 
through simulations. 

The QRA process increases the RMS noise per pixel at 
the level ~ 1% for a/q~2. The distribution of the RMS 
induced by the QRA process depends on the scanning law 
and resembles the distribution of white noise in the map. 
At low and middle ecliptical latitudes, for Planck/LFI 
the RMS decreases toward the ecliptical poles scaling ap- 
proximately as vsin^ where 9 is the ecliptical colatitude. 
The details of the distribution of the QRA additional noise 
on the map depend on the frequency channel, the location 
of the horn inside the focal plane and the scanning law. We 
have derived the parameters that characterize this depen- 
dence for Planck/LFI. We find that the skewness and the 
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N s 


8n s 


(J iV B 




1 


0.3414 


0.2887 


1.1825 


2 


0.2183 


0.2041 


1.0696 


3 


0.1726 


0.1667 


1.0359 


1 


0.1482 


0.1443 


1.0269 


5 


0.1321 


0.1291 


1.0230 


6 


0.1205 


0.1179 


1.0227 


10 


0.0930 


0.0913 


1.0187 


15 


0.0753 


0.0745 


1.0108 


30 


0.0529 


0.0527 


1.0033 


60 


0.0373 


0.0373 


1.0017 



Table A.l. Table of values of 5n s 
function of AL. 



<7jv a and o>n s as a 



kurtosis measured all over the map are modified up to a 
few percent for a/q ~ 2 when the 1// Q noise is subtracted 
with dcstriping algorithms. Although not large, this effect 
is not negligible compared to the sensitivity of the kur- 
tosis index recovery for a pure CMB sky. On the other 
hand, it is much less relevant than residual foreground 
contamination which accurate modelling turns to set the 
fundamental uncertainty on the recovery CMB anisotropy 
high order moments. 

Finally, we find that the quantization introduces a 
power excess, Cf x , that, although related to the instru- 
ment and mission parameters, is weakly dependent on the 
multipole £ at middle and large £ and can be quite accu- 
rately subtracted. For a/q ~ 2, the residual uncertainty, 
ACf x , implied by this subtraction is of only ~l-2% of the 
RMS uncertainty, AC™ mse , on C\ ky reconstruction due to 
the noise power, C™ OTSe . 

Only for £<30 the quantization removal is less accu- 
rate; in fact, the 1// noise features, although efficiently 
removed, increase C£ mse , AC £ " oise , Cf and then ACf* . 
Anyway, at low multipoles Cf y > AC^ noisc > ACf x and 
the uncertainty introduced by the QRA effect is therefore 
in any case much less than the unavoidable uncertainty 
due the cosmic variance. 
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Appendix A: calculation of ajs s 

The ol ty s coefficient in eq. (|13|l is defined as the ratio be- 
tween the standard error for the QRA process after N s 
averages Sn s and the RMS for the same process: 



(A.l) 



where (Tn s is the RMS of the QRA error distribution after 
N s averages: UN a (X), (5jv s is obtained from the definition 
of standard error: 



Ill dtuN ^ 

V 2tt j 1 



0.6827. 



(A.2) 



We are interested in determine the central moments and 
the integral of un s (5n b ) i- e - to determine the characteristic 
function of ujv 3 (<5at s ), from the characteristic function for 
Sjv a = NgS^. Where, Sjy a represents the distribution of 
the sum instead of the average of N s random variables. 
Taking in account that both distributions are even and 
for N s = 1 they coincide with the top-hat distribution: 
u(x) = 1, for \x\ < 1/2, or for |x| > 1/2, it is possible to 
demonstrate that the characteristic function of UN a {SN a ) 
is: 



sync 



2N, 



(A.3) 



where sync(x) = sin(cc)/a;. From eq. (|A.3|I it is imme- 
diate to recover the central moments of un s (5n s )- In 
particular it is relevant to note that all the odd mo- 



ments are null, while: mi = o jy / q ■ 
{5N S - 2)/(24(W s 3 ) = kurtosis/g 4 . 



1/(12AQ and 



m 4 



The evaluation of Sn s requires to solve eq. 
is obtained by the following formula: 

In s {$nJ = R(uj max )+ 

4(5 



H. This 



^ JIT""* s y nc (<"*«-.) [ s y nc (5^) 



(A.4) 



du) , 



where w max represents a cut-off in an otherwise indefinitely 
extended integration range. The residual i?(a; max ) is up- 
per bounded by: i?(w max ) < 4(2N S ) N °/ (ttN s lo^iJ) which 
allows to determine w max 
integration e is given: 



once the absolute accuracy of 



= 2AL 



"I l/JVs 



7reAT R 



(A.5) 



Note that it would be also possible to solve eq. I|A.2(1 
through direct integration of Ujv B (<5jv 3 )j but this is not a 
satisfactory method since the integral converges slowly for 
N a > 10. 

Values of evjv a as a function of N s from equa- 
tions IjA.ip and (|A.2|I are tabulated, with four digits of 
accuracy and for N s = 1, 2, 3, 4, 5, 6, 10, 15, 30, 60, in 
the last column of Tab. IA.l1 While the values of ^jv a , and 
a tv s are tabulated in the second and the third columns of 
the table respectively. For N s > 30 the following relation 
holds: 

log 10 (ajv 3 - 1) « -0.923359 • log 10 N s - 1.129003. (A.6) 

with an accuracy compatible with the accuracy of the tab- 
ulated values. 
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Appendix B: the ADC quantization problem 

In this paper we have extensively discussed the effect of 
the software quantization of the data before lossless com- 
pression. This quantization is only the last step for the 
on-board data-processing which may affect the data qual- 
ity, being the data acquired by the on-board computer by 
an analog to digital converter (ADC) which introduces at 
least three sources of noise: i) the intrinsic ADC quanti- 
zation noise proportional to the ADC quantization step 
9adc; H) the read-out noise oron; iii) non linearities due 
to residual differences in (?adc as a function of the input 
value seen by the ADC. Detailed analyses of these prob- 
lems are outside the scope of this paper and will be ad- 
dressed in forthcoming works. However, in the context of 
this study on the signal digitisation effect it is relevant to 
discuss the conditions under which the ADC quantization 
effect is not critical compared to the software quantization 
effect. 

The impact of the ADC quantization depends on the 
steps of the on-board data processing. Assuming (a) to ac- 
quire separately the sky and the reference signal with an 
ideal ADC (i.e. considering negliglible the sources of noise 
i) and ii)), (b) to have a perfectly balanced radiometer (i.e. 
both the sky and load detectors have the same weight), (c) 
the software quantization, and then the total QRA error, 
to be independent from the ADC quantization, and (d) to 
have ATh ADC samples that are averaged to obtain a given 
instrumental sample (i.e. the hardware sampling rate is Ah 
times higher than the software sampling rate), then the 
global RMS quantization error, a g , including both hard- 
ware and software quantization, is given by: 



Ah 



<7ADC 



(B.l) 



where a q is the RMS of the software quantization error. 
This relation is substantially correct at the level of TODs 
as well as at the level of the maps, since the repeated 
observation of the sky pixels and the averaging determined 
by the scanning strategy operate in the same way in the 
case of data affected by software quantitazion alone or by 
hardware plus software quantization, analogously to the 
considerations of Sect. 2.1. Clearly, a g is larger than the 
expected QRA error reported by eq. which represents 
a very good approximation in the limit q 3> qxac and/or 
for large values of Ah- 

In order to understand the conditions under which the 
ADC quantization effect is not critical compared to the 
software quantization effect we focus on the impact on 
the angular power spectrum estimation. At middle and 
large multipoles the power spectrum of the global and the 
software quantization errors are proportional to the square 
of the corresponding RMS quantization error, i.e.: 



£tQRA,g 



QRA 



(B.2) 



By following an approach analogous the that described 
in the second part of Sect. 4, we require that 



C QRA, B _ pQRA ^ 1 RMS^QRAj 



QRA 



a a 



QRA 



(B.3) 



Differently from the case discussed in Sect. 4, we have 
now to search for a law, L/>, which gives at each £ the 
minimum value among the ratios RMS[Cfgo]/Cfgo found 
for the set of simulations considered here. We find that 
the approximate law 



0.46f 



-0.5 



(B.4) 



is a quite good approximation for the whole range of mul- 
tipoles relevant here. 

By using the above relation between a g and o~ q and this 
expression for Li, the above condition is satisfied provided 
that 



J_ fgAPoy , 0.23l°- 5 
N h V q J ~ A 



(B.5) 



In the case of Planck/LFI, Ah ~ 50 and even in the 
worst case of q ~ 0.5 mK and ^adc — 1.2 mK, this con- 
dition is satisfied at £^30 with a value of A ~ 20, which 
corresponds to a deviation from the value of expectation 
of the angular power spectrum of pure software quanti- 
zation error of only few per cent of its RMS. Of course, 
given an accurate description of the ADC quantization we 
can easily include it in the data analysis. Anyway, this 
estimate shows that its effect is typically very small. 

Finally, it is worth to note that just replacing 2/ Ah 
with 1 /Ah in eq. IjB.lfl accounts for a different conceptual 
acquisition scheme where the sky and reference signals are 
differenced before to perform the ADC quantization. This 
scheme has a better propagation of the quantization error 
than the scheme assumed in eq. IjB.ljl . implying that the 
above condition can be satisfied with a value of A two 
times larger. 

Appendix C: Statistical moments of simulated sky 
maps 

We report here a comprehensive tabulation of the statis- 
tical moments of the CMB simulated map adopted in this 
work compared with the statistical moments of simulated 
maps of the most relevant Galactic and extragalactic fore- 
grounds. Our maps of Galactic foregrounds are simulated 
according to Maino et al. (2002), that of extragalactic 
source fluctuations according to the model by Toffolatti 
et al. (1998), while for thermal Sunyaev-Zeldovich effects 
(SZ) from clusters of galaxies we have used the template 
available at the MPI web site 111. It is infact interesting 
to compare the modification on the statistical moments of 
a pure CMB anisotropy sky due to the digisation effect 
with that due to the foreground contamination. 
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The HEALPix scheme (Gorski et al. 1998) is here 
adopted and different resolutions, identified by the pa- 
rameter -/Vgidc (12 x -ZV s 2 ido is the number of pixels in the 
map), are considered. 

We consider here maps with and without beam convo- 
lution (FWHM of 33' as appropriate to the LFI 30 GHz 
channels). The input maps have been first simulated at 
N s ide = 1024 without beam convolution. Convolution and, 
possibly, degradation are then applied. 

We report the statistical moments referring separately 
to each component and to combination of CMB and fore- 
ground maps, by including or not masks to exclude regions 
at low Galactic latitudes and/or high signal pixels, as in- 
dicated in the tables. In fact, the regions at low Galactic 
latitudes have to be avoided in CMB anisotropy analysis 
as well as pixels significantly contaminated by SZ effects 
or extragalactic sources have to be previously detected 
and removed. More precisely, for uniformity reasons, we 
apply the same 5cr-clipping threshold, a being the rms 
of CMB fluctuations at the considered resolution for the 
corresponding case (including or not beam convolution), 
independently of the considered kind of map composition. 

Note that, while in the cases of unconvolved maps with 
large skewness and kurtosis indices the convolution de- 
creases both these estimators, in the cases of unconvolved 
maps with moderate or small skewness and kurtosis in- 
dices the convolution may produce a weak increase of 
these estimators because the higher order moments are 
relatively less reduced than the variance. 
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